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The dynamics of an open quantum system is usually studied by performing a weak-coupling and 
weak-correlation expansion in the system-bath interaction. For systems exhibiting strong couplings 
and highly non-Markovian behavior this approach is not justified. We apply a recently proposed 
correlated projection superoperator technique to the model of a central spin coupled to a spin bath 
via full Heisenberg interaction. Analytical solutions to both the Nakajima-Zwanzig and the time- 
convolutionless master equation are determined and compared with the results of the exact solution. 
The correlated projection operator technique significantly improves the standard methods and can 
be applied to many physical problems such as the hyperfine interaction in a quantum dot. 

PACS numbers: 03.65.Yz, 42.50.Lc, 03.65.Ta, 73. 21. La 



I. INTRODUCTION 

Some of the key features in the dynamics of open quan- 
tum systems [l| are phenomena such as relaxation, deco- 
herence and the buildup of correlations and entanglement 
due to the interaction of the open system with its envi- 
ronment. This behavior is generic in the sense that it is 
observable for many different kinds of environments and 
microscopic interactions. For many physical systems it is 
justified to assume that the coupling to the environment 
is weak (Born approximation) or that correlations in the 
bath decay quickly with respect to the typical timescale 
of the system's dynamics (Markov approximation). In 
this case it is often possible to construct a generator for 
the dynamics which is in Lindblad form 0,9- 

In general, however, memory effects in the bath can- 
not be neglected and the Markov assumption is not ap- 
plicable any more. Such systems are said to exhibit non- 
Markovian behavior. This can be due to strong system- 
environment couplings L correlations and entangle- 
ment in the initial state @, 0] , finite reservoirs d, [1| , or 
due to coupling to environments at low temperatures or 
to spin baths [13, [HI [ill ■ Also, heat transport in nano- 
structures [l3| has been shown to exhibit non-Markovian 
behavior. 

A powerful tool for dealing with such systems is pro- 
vided by the projection operator techniques [l^ . ITsj 
which have been introduced by Nakajima [l6[, Zwanzig 
and Mori [l^. These techniques are based on the 
introduction of a projection superoperator V which acts 
on the density operator p of the total system and projects 
onto the so-called relevant part Vp. The latter represents 
a certain approximation of p, and this procedure leads to 
a simplified effective description of the dynamics through 
a reduced set of relevant variables. Stated differently, the 
projection V expresses the elimination of certain degrees 
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of freedom which are considered negligible for the treat- 
ment of the open system's dynamics. 

Here, we will consider two different approaches which 
both lead to a closed equation for the dynamics of the 
relevant part of p. The first one leads to the Nakajima- 
Zwanzig equation [l^, [13] , an integrodifferential master 
equation for Vp containing a memory kernel. The second 
approach eliminates the integration over the system's his- 
tory and leads to the time-convolutionless master equa- 
tion which is local in time, but involves a time-dependent 
generator [HIIO, HH [H [M HI ■ In both cases the result- 
ing master equation is used as a starting point for a sys- 
tematic expansion in powers of the system-environment 
interaction. 

Usually the projector is chosen such that Vp = 
{trsp) po — Ps ® Pa, where pa is some fixed environ- 
mental state, e.g. a thermal equilibrium state, and where 
tr^; denotes the partial trace over the environment. The 
superoperator V projects the total state p onto a tensor- 
product state, i.e., onto a state without any statistical 
correlations between the open system and its environ- 
ment. If correlations remain small, this ansatz is justi- 
fied and the perturbation expansion is usually carried out 
only to leading order in the perturbation (Born approx- 
imation). However, for many physical systems exhibit- 
ing strong correlations, the Born approximation might be 
inadequate such that one also has to take into account 
higher orders of the expansion. The calculation of higher 
orders is limited by the increase of mathematical com- 
plexity, though, and furthermore, the expansion may not 
converge uniformly in time, such that higher orders may 
diverge on longer timescales ^25'] . 

A different strategy for taking strong correlations into 
account is to introduce a superoperator V that projects 
onto a correlated system-environment state, i.e., onto 
a state that contains certain statistical correlations be- 
tween the open system and its environment [1^, [13, [H, 
[H, [13, [m [111 ■ We will refer to such P as a correlated 
projection superoperator. In the present article we will 
discuss in detail the application of the correlated projec- 
tion operator method to a two level system interacting 
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with a spin bath. This system is particularly interesting 
since it models the hyperfine interaction of an electron 
confined to a quantum dot with its surrounding nuclei 

We organize this paper as follows. In Sec. |TT] we re- 
view the standard projection operator technique and 
formulate the Nakajima-Zwanzig as well as the time- 
convolutionless master equation. Also, we briefly recapit- 
ulate how to extend the technique to correlated projec- 
tion superoperators. In Sec. IIIII we apply this method to 
a simple spin-bath model and construct projectors on the 
basis of symmetry considerations that vastly improve the 
performance of the projection operator technique. We 
compare the resulting solutions of the Nakajima-Zwanzig 
and the time-convolutionless master equation in second 
order to an exact solution of the model. Discussion and 
conclusions are then provided in Sec. IIVI 



II. PROJECTION OPERATOR TECHNIQUES 

In general, one considers an open system S which is 
coupled to some environment E. The accordant Hilbert 
spaces are denotes by Tig and He, respectively, and the 
state space of the total system is given by the product 
Ti = Hs ® He- The Hamiltonian of the total system is 
denoted hy H ^ Hq + Hi, where Hq is the unperturbed 
part of the Hamiltonian and Hi represents the interaction 
between system and bath. A general state of the total 
system is given by a density matrix p. The partial traces 
over the system S and the environment E are denoted by 
tr5 and tr^;, respectively. The reduced density matrix of 
the open system is thus given by ps = tr^p. 



A. The standard projection superoperator 



B. Nakajima-Zwanzig equation 

Starting from the von Neumann equation in the inter- 
action picture with respect to Hq, 



dt 



p{t)^-i[Hi{t),p{t)]=C{t)p{t), 



(2) 



one can derive a closed equation for the projection Vpit) 
by inserting V and Q in front of p on both sides. 



dt 



Vpit) 



dti ICit,ti)Pp{ti). 



(3) 



The superoperator 



)C{t,ti) = VC{t) Texp 



dt2QC{t2) 



QC{ti)V (4) 



is referred to as the memory kernel or the self-energy, 
and the operator T denotes chronological time-ordering. 
Eq. ([3]) is called the Nakajima-Zwanzig (NZ) equation 
and describes non-Markovian behavior of the total sys- 
tem through the memory kernel ((4]). In general, there is 
an additional term proportional to Qp{Q) on the right- 
hand side of the NZ equation which we have omitted 
here for simplicity, supposing an initial state that satis- 
fies 'Pp(O) — p{0). Since IC{t, ti) is usually a very compli- 
cated operator, it is customary to perform a perturbation 
expansion in powers of Hi. Under the condition 



rc{t)r = 



(5) 



the lowest-order contribution is given by the second or- 
der, 



dtirc{t)c{ti)Pp{ti 



(6) 



and higher orders are obtained by expanding the time- 
ordered exponential of the memory kernel (|4]). 



The starting point of the projection operator tech- 
niques is the introduction of a superoperator "P which 
acts on the total system's density matrix, and which is 
defined by 

Vp= (tlEp)® PO, (1) 

where po is some fixed state of the environment. The map 
V satisfies the condition of a projector, namely — V, 
and is therefore referred to as a projection superopera- 
tor, being a map acting on operators. The complemen- 
tary map is defined via Q = I — V , where / denotes the 
identity. Note that Vp contains all information about the 
open system in the sense that for the expectation value 
of any observable Os of the open system the relation 
tr{Osp} = tT:{OsVp] holds. 



C. Time-convolutionless master equation 

An alternative way of deriving an exact master equa- 
tion for the relevant part of p is to remove the dependence 
of the system's dynamics on the full history of the system 
and to formulate a time-local equation of motion, which 
is given by 

lvp{t)=JC{t)'Ppit). (7) 

This equation is called the time-convolutionless (TCL) 
master equation, and IC{t) is a time-dependent superop- 
erator, which is referred to as the TCL generator. As for 
the NZ equation, in general there is also a term propor- 
tional to Qp{0) on the right-hand side of Eq. ([7]), which 
vanishes provided one chooses a factorized initial state. 



3 



Like for the NZ equation one can carry out a pertur- 
bation expansion of the TCL generator in powers of Hj. 
The various orders of this expansion can be expressed 
through the ordered cumulants [H, [s^, [1^ of the 
Liouvillc superoperator C{t). The second-order contri- 
bution reads 

/C2(<) = / dhVC{t)C{ti)P, (8) 
Jo 

such that the second-order TCL master equation takes 
the form 

j^Vpit) = ^ dhVC{t)C{h)Vp{t), (9) 

which should be contrasted to the second-order NZ equa- 
tion One can formulate a simple set of rules that 
enable one to write down immediately these expressions 
and the corresponding expressions for all higher orders 
The TCL technique has been applied to many phys- 
ical systems. Examples, which include the determina- 
tion of higher orders of the expansion, are the damped 
Jaynes-Cumings model (3 Til, q uantum Brownian motion 
and the spin-boson model [38j, the spin-star model psj . 
relaxation processes in structured reservoirs [27j , and the 
dynamics of the atom laser (s^ . Applications of the TCL 
method to systems relevant for quantum information pro- 
cessing may be found in [13, Sli] . Recently, the perfor- 
mance of the TCL approach has been studied for the 
model of a qubit coupled to a spin bath via Ising inter- 
action [i^l, and for a general open system coupled to a 
bath through pure phase noise [4^. 

We emphasize at this point that the TCL master equa- 
tion describes non-Markovian dynamics although it is lo- 
cal in time. All memory effects are taken into account 
by the explicit time-dependence of the TCL generator. 
It should also be pointed out that both the NZ equation 
^ and the TCL equation ([7]) are exact in the sense that 
they are equivalent to the equation of motion gained from 
the full system Hamiltonian. 

It is important to realize that the NZ and the TCL 
techniques lead to equations of motion with entirely 
different structures. Therefore, also the mathematical 
structure of their solutions are quite different in any given 
order j!!]. It turns out that in many cases the degree of 
accuracy obtained by both methods is of the same or- 
der of magnitude. In these cases the TCL approach is 
of course to be preferred because it is technically much 
simpler to deal with. 

D. Correlated projection superoperators 

The performance of the projection operator techniques 
depends on the properties of the microscopic model under 
study, in particular on the structure of the environmental 
correlation functions. However, it also depends strongly 
on the choice of the superoperator V . Several extensions 



of the standard projection JT]) and modifications of the 
expansion technique have been proposed in the literature 
(see, e.g., (H,!!!!!). 

A characteristic feature of the projection defined by 
Eq. (U]) is given by the fact that it projects any state 
p onto a tensor product ps ® po that describes a state 
without statistical correlations between the system and 
its environment. This is not the only possible choice. In 
fact, a general class of projection superoperators can be 
represented as follows, 

Vp^^trE{A,p}(g)B„ (10) 

i 

where {Ai} and {Bi} are two sets of linear independent 
Hermitian operators on He satisfying the relations 

trE{B,A,} = (11) 
^(tr^BOA, = Ie, (12) 

i 

J2Af<»B, > 0. (13) 

i 

The map defined by Eq. (jlOp projects in general onto cor- 
related (non-factorizing) system-environment states and 
may thus be referred to as a correlated projection super- 
operator. It can be shown [2^ that Eq. (fTO)l represents 
the most general form of such a correlated projection 
superoperator under certain natural physical conditions. 
These conditions demand that "P is a completely posi- 
tive and trace preserving map (quantum channel) that 
operates on the environmental variables. Equation pip 
guarantees that 7^ is a projection superoperator, Eq. (fT^ 
ensures that V preserves the trace (normalization), while 
Eq. (fT3)) is equivalent to the condition of complete posi- 
tivity (T denotes the transposition) . We remark that the 
class of projections defined by Eq. pn]) yields a natural 
generalization of the Lindblad equation to the reg ime of 
strong non-Markovian quantum dynamics [26l l28j 

A typical example for a correlated projection that is 
relevant for the present paper is obtained by the choice 

A = n„ = -^n„ iv, = tr£;n„ (w) 

where the 11^ are ordinary projection operators on TCe 
which form an orthogonal decomposition of the unit op- 
erator: 

n,n, = Jy-Hj, n| = n„ ^ n, = /b. (15) 

i 

One easily checks that with this choice all conditions for- 
mulated above arc indeed satisfied. Further examples are 
discussed in Refs. [26, 27, 28, 31, 32]. 

As mentioned already in Sec. Ill C[ a homogeneous NZ 
or TCL equation of motion presupposes a tensor product 
initial state if one uses the standard projection superop- 
erator ([T]). However, this is no longer true for the corre- 
lated projection defined by Eq. (fTU)) . In fact, the general 
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condition for the absence of an inhomogeneous term in 
the NZ equation ^ or the TCL equation ([7]) is given 
by Vp{0) — p{0). According to Eq. (fTO|l this condition 
does not require that p(0) be a factorized state. Hence, a 
great advantage of the correlated projection superopera- 
tors is given by the fact that they allow the treatment of 
correlated initial states by means of a homogeneous NZ 
or TCL equation 26, 4^. 

Once V is chosen, the dynamics of the open system is 
uniquely determined by the dynamical variables 

p,{t)^tlE{Ap{t)}. (16) 

The connection to the reduced density is simply given by 

i 

and the normalization condition reads 

trsPs{t)=Y,tl■sp^{t)^l. (18) 

i 

The reduced density matrix is hence uniquely determined 
by a set of (unnormalized) operators pi{t). With the 
help of the NZ equation ([3]) one finds that the dynamics 
of these operators is governed by a coupled system of 
integrodifferential equations, 

d /■* 

J^P^{t) = "^^1 IC,j{t,h)pj{h), (19) 

where the superoperators Kijlt,ti) are defined through 
their action on an arbitrary system operator Os, 

/C„ (i, h)Os = trs {A,lC{tM){Os ® Bj)} . (20) 

On the other hand, the TCL equation d?]) leads to a cou- 
pled system of time-local differential equations, 

±p,{t)^Y.^,,{t)p,{t), (21) 

3 

with superoperators defined by 

lC^,{t)Os = trs {A,K{t){Os ® B^)] . (22) 

A crucial step for a successful applications of the cor- 
related projection operator technique is the construction 
of an appropriate projection superoperator V. A use- 
ful strategy for this construction is to take into account 
the known conserved quantities of the model under study. 
Let us suppose that C is a such conserved quantity. Then 
a good choice for 7-' is a projection that leaves invariant 
the expectation values of C, i.e., tr{Cp} = tr{CVp}. Re- 
quiring this to hold for all states p we get the condition 

V^C = C, (23) 

where denotes the adjoint of V with respect to the 
Hilbert-Schmidt scalar product [48|. This equation rep- 
resents a condition for the projection superoperator V 



on the basis of a known conserved quantity of the under- 
lying model. It ensures that the projection superopera- 
tor leaves invariant this quantity and that the effective 
description respects the corresponding conservation law 

mm- 

III. APPLICATION TO A SPIN BATH 

A. Description of the model 

We want to apply the correlated projection operator 
formalism to a spin-bath model, which is defined by the 
Hamiltonian 

N 

^= (24) 
fc=i 

The central spin, representing the open system, has an 
energy splitting ojq and is coupled to the fc-th bath spin 
via the coupling constant A^. Here, <t and cr^*^^ are the 
Pauli operators of the central spin and of the k-ih bath 
spin, respectively. 




FIG. 1; (Color online) Illustration of the model with Hamil- 
tonian (|24[) . A central spin cr interacts with a bath of N spins 
cr'*' through the coupling constants Ak- 
in the following we will restrict ourselves to initial 
states of the form 

p{0) = ps{0)®^lE, (25) 

where Ie is the unit matrix in the state space of the spin 
bath. Initially the bath is thus in a completely unpo- 
larized (infinite temperature) state. This choice is made 
here only for simplicity; the present method also allows 
the treatment of polarized bath states and of correlated 
initial states. 

Our aim is to illustrate the TCL and the NZ expansion 
for a correlated projection superoperator. To compare 
the results with the exact solution of the full von Neu- 
mann equation corresponding to the Hamiltonian (j24p 



5 



we take the coupling constants to be independent of k, 
i.e., = A. The model can then be solved analytically 
as follows. First one notes that the 3-coniponent of the 
angular momentum of the total system, 



= 77^3 + -h, 



(26) 



is conserved under the time-evolution given by the Hamil- 
tonian p4)) . Here, J3 denotes the 3-component of the 
angular momentum of the bath which is defined by 



(k) 



(27) 



For Ak = A, also the square of the bath angu- 
lar momentum is conserved. We introduce basis states 
I j, m) which are defined as simultaneous eigenstates of 
and J3 [4^. If N is even, j takes on the values 
j = 0, 1, 2, . . . , f , while j = i, f , . . . , f for odd N. For 
a given value of j, the quantum number m takes on the 
values m = -j, -j + 1, . . . , +j. 

With the help of the conserved quantities one easily 
solves the time-dependent Schrodinger equation, employ- 
ing the fact that the two-dimensional subspaces spanned 
by the states |-|-) (8) \j,m) and |— ) (g) \j,m + 1) are in- 
variant under the time-evolution. Here, the states |±) 
denote the eigenstates of the central spin, (T3|±) = ±|±). 
We further introduce the quantities 



where 



N 



N 



N 



(28) 



(29) 



represents the number of times the angular momentum j 
appears in the decomposition of the bath Hilbcrt space 
into irreducible subspaces of the rotation group. Hence, 
p{j) may be interpreted as the probability of finding a 
certain angular momentum j in the initial state of the 
bath. Finally, we define the frequencies 

n±{m) ^ ±wo + 4A(±TO+ 1/2), 
t^±{j,m) = \ + 4A2[j(j + 1) - m(m ± 1)]. 



With the help of these definitions the upper state popu- 
lation of the central spin can be written as 



P+{t) = {+\ps{t)\+) 



E 



2j + l 



cos'^ifJ' +{j,m)t] 



4Ai+(j,"i) 



sin^[^+(j, m)t] , (30) 



where we have used for simplicity the initial condition 
P+(0) = 1. The lower state population is found from 



p_^t) = {-\ps{t)\-) = 1 - P+{t). The coherence of the 
central spin is determined by the formula 



p+^it) = {+\psit)H 



pU)e 



iujQt 



jm 



23 + 1 



r / N 1 iO.+ (m) 
cos[p+{3,m)t\ - sin[p+{j, m)t\ 

cos[p.-{j, m)t] + sin[/i_(j, m)t] 

2p.^{j,m) 



B. Correlated projection operators 



According to Sec. IH Dl a good candidate for a corre- 
lated projection supcroperator is a map that does not 
change the expectation values of the known conserved 
quantities. We start by using a projection supcropera- 
tor V that leaves invariant the 3-component of the total 
angular momentum, i.e., that satisfies [see Eq. ([23| ] 



^tjtot 



rtot 
-^3 7 



(31) 



We define bath operators n,„ which project onto the sub- 
space spanned by the eigenstates of J3 corresponding to 
the eigenvalue m = — appropriate pro- 
jection operator with the property (j3ip is then given by 



7'p = ^tr£;{H„p}( 



where the number 



1 



-H. 



^Pni'Si 



1 



-Tlr 



N 



(32) 



(33) 



represents the degree of degeneracy of the eigenvalue m 
of J3. One can easily check that this projection operator 
satisfies all conditions formulated in Sec. IH Dl 

The state of the central spin is thus determined by the 
set of densities Pm(t)- The projection Vp represents a 
correlated system-environment state in which each pm is 
correlated with the bath state Hm/iVm which describes a 
state of maximal entropy under the constraint of a given 
value m for the 3-component of the angular momentum. 
According to Eq. (fTT]) the reduced density matrix is ob- 
tained from 



PS 



(i) = ^p„(t). 



and its normalization condition reads 

trs Ps{t) = ^^rs Pm{t) = 1. 



(34) 



(35) 



Finally, the matrix elements of the pm{i) are defined by 

Pi"(t) = {±\pm{t)\±), (36) 



(37) 
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It should be emphasized that the above procedure is also 
applicable if the coupling constants A^. are not equal to 
each other, since J^°^ is then still a conserved quantity. 
We remark that the standard procedure with a projection 
operator of the form of Eq. ([T]) with pQ = Ie would 
give jtot ^^cr^. 



Interaction Hamiltonian 



We write the Hamiltonian (|24p in the form H = Hq 
Hi, where 



UJn 



is regarded as the unperturbed part and 
Hi = 2 A (f7+ J_ +cr_J+) 



(38) 



(39) 



as the interaction Hamiltonian. Note that the unper- 
turbed Hamiltonian Hq is not a sum of free Hamiltonian 
operators for system and environment, and that the per- 
turbation Hamiltonian contains only the flip-flop part of 
the interaction, which describes the flip of the central 
spin accompanied by a back-flip of one of the bath spins. 

We transform to the interaction picture with respect 
to Hq which leads to the interaction Hamiltonian 



Hi{t) = 2Ae'^ 



ujo-2A)t 



-t-H.c. 



We will always represent our results in the interaction 
picture with respect to the Hamiltonian ^cfs, i.e. in the 
rotating frame of the central spin. The results in the in- 
teraction picture relative to Hq are transformed back to 
the rotating frame of the central spin through the pre- 
scription 



Pra{t) 



(40) 



Obviously, this back-transformation leaves the popula- 
tions P±{t) unchanged, while the coherences pick up a 
phase factor. 



—AiAmt m 



(t). 



2. Master equation 



For the projection defined by Eq. (15^ the NZ master 
equation (fTO|) takes the following form in second order. 



d_ 

di 



"'0 



X tris <^ H„ 



Hi{t), 



Hliti),pm'iti) ® 



N„ 



-H„ 



The right-hand side of this equation is easily calculated 
to yield [13] 



dt 



Pm{t) = f dh 

Jo 



{N \ 

a+p,n+l{tl)(T-y— + 771+ lj2cOs[17+(77l)(t - tl)] 

/N \ 
-a^Pm~i{ti)a+y— - 777 + lj2cos[r2„(777)(t - tl)] 



N 



-(7+fT_p,„(il)(^ — 



N 



-a-a+p„i{ti)(^ — 



me 



^iSl_ (m)(t-ti) 



N 



fN 



+ 777 e 



-iSl_ {m)(t-ti) 



(42) 



One easily checks that this master equation preserves the 
Hermiticity and the trace of the reduced density matrix 
ps{t) which is determined by Eq. (p4|) . Note that the dy- 
namics couple the densities Pm{t) to Prn±i{t)- The corre- 
sponding second-order TCL master equation Eq. (|2ip is 
obtained from Eq. (|42)l by just replacing Pm{ti) by Pm{t) 
on the right-hand side. 



3. Coherences 

The coherence of the central spin is given by the sum 
of the coherences of the , 



From the initial condition ((25l) we find that 

2^ 



-(0) 



Jp+-io). 



(43) 



(44) 



The dynamics of the coherences p'^_ (t) is determined by 
the master equation (H^ which leads to 



B+(777)e'"+('")(*-*i) 

PT-(ti),(45) 



dt'--^'^ - -Jq 



(41) where 



+B-(77i)e"*^^-(™^(*~*i) 

B±(777) = 4A^ 



N 



+ 777 



(46) 



Note that the dynamics does not couple the different 
components p^_{t). The integrodifferential equation 
can now be solved with the help of a Laplace trans- 
formation. The corresponding TCL equation is obtained 
from Eq. P5)l by replacing p™_(ii) by p+_(i) on the 
right-hand side. Solving the resulting equation yields 



p+_(t) = p+_(0) ^ exp [-AtAmt - A^^i)] , 



(47) 
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where 



B+{m) 
B-{m) 



-ir2_ {m)t 



B+{m) B^{m) 



ri+(rn) r2_(rn) 



(48) 



and the TCL equation reproduce the exact solution very 
well in the perturbative regime. Moreover, both are ca- 
pable of describing the position, the height and the width 
of the partial revivals of the coherence, which appear be- 
cause of the finite size of the bath and the associated 
finite recurrence times. For large integration times, only 
the positions and the width of the revivals are repro- 
duced correctly (see Fig. H]). We note that the NZ and 
the TCL method yield nearly identical results; they are 
hardly distinguishable on the scale of the shown figures. 
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FIG. 2: (Color online) Coherence of the central spin. Con- 
tinuous line: Exact solution. Dashed line: Second-order NZ 
approximation. Parameters: A'' = 101 and a — 0.1. 



FIG. 4: (Color online) Coherence of the central spin. Con- 
tinuous line: Exact solution. Dashed line: Second-order TCL 
approximation. Parameters: A'^ — 101 and a = 0.1. 
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FIG. 3: (Color online) Coherence of the central spin. Con- 
tinuous line: Exact solution. Dashed line: Second-order TCL 
approximation. Parameters: N — 101 and a = 0.1. 



The comparison with the exact solution shows that 
both methods correspond essentially to an expansion in 
the dimcnsionless quantity 



a 



2AN 

Wo 



(49) 



We compare the results of the NZ and the TCL method 
with the exact solution in Figs. [Hand [31 where the coher- 
ence is normalized such that (0) = 1. Both the NZ 



4. Populations 
The populations are given by 

p±(i)-E^± 

m 

The initial condition (f25|) now yields 

^r(0) = ^' P!"(0) = 0. 

The master equation leads to 

d 
dt 

' 'N 



(50) 



(51) 



— +m + l] cos[n+{m){t - ii)]P™+i(ii) 



^-m)cos[17+(m)(i-<i)]^^+(ii) 



(52) 



and 



d_ 
di 



P!^{t) = 8A^ [ dh 
Jo 

N \ 

y + mj cos[l}_(m)(t-<i)]P™(ti) 



(53) 
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Note that, in contrast to the coherences, this equation 
does couple different components P™. However, the so- 
lution can still be constructed by noting that 



-[P™(t) + P™+i(t)] =0. 



(54) 



With the help of this equation one shows that the to- 
tal angular momentum Jg"* is exactly conserved under 
the dynamics given by the master equation. Namely, we 
deduce from Eq. ([54)l that 



[P^{t)+P"'+\t)] =0. (55) 



-tr{jftpp(t)} 



-t(- 



X 10"' 




X 10" 
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FIG. 5: (Color online) Population of the central spin ac- 
cording to the exact solution (continuous line), and to the 
second-order TCL approximations using the correlated pro- 
jection (broken line) and the standard product-state projec- 
tion (dotted line). Parameters: A'^ = 101 and a = 0.5. 



FIG. 6: (Color online) Population of the central spin accord- 
ing to the second-order TCL approximation (top) and to the 
exact solution (bottom). Parameters: = 101 and a = 0.5. 



their height is not (see Fig. [6]). We find again that the 
NZ result is hardly distinguishable from the correspond- 
ing TCL result, and hence refrain from showing it in a 
separate figure. 

The standard procedure that uses the product-state 
projection operator ([1]) with po 2"^/^; yields the upper 
state population 



1 -I- exp 



(1 — COSWo^) 



(58) 



which is also plotted in Fig. [5] We see that the standard 
projection operator technique leads to a rather bad ap- 
proximation which is neither able to represent correctly 
the damping nor the revivals of the populations. 



Equation ((54)) is used to eliminate P™+^ from the 
right-hand side of Eq. ([52]) to get a closed equation for 
P™ which can again be solved with the help of a Laplace 
transformation. The corresponding TCL equations are 
solved in a similar manner yielding 



where 



2N 



Y + m + 1 ^ Y !I^;g-A^°p(t) 



iV+ 1 



iV+ 1 



SA^iN +1) 

nlim) 



(1 - cos[n+{m)t]) . 



(56) 



(57) 



As demonstrated in Fig. \5\ this result gives an excellent 
approximation of the exact dynamics in the perturbative 
regime. For longer integration times, the position and 
the width of the revivals are reproduced very well, while 



Complete projection onto the angular 
momentum manifolds 



For the present model we can construct another corre- 
lated projection superoperator which takes into account 
both conserved quantities, namely Jg"* and J^. This 



projection is given by 



jm 



trE{TljrnP} 



-n. 



-n, 



(59) 

where Nj has already been defined in Eq. . The Hjm 
denote the ordinary projection operators that project 
onto the common eigenspaces of J3 and with cor- 
responding eigenvalues m and j(j + 1). The map ([59]) 



obviously leaves invariant Jg"* and and fulfills again 
all conditions formulated in Sec. IIIDI 

The second-order NZ master equation corresponding 
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to the projection ([59]) reads 



dt 



p,^{t)^AA' / dti 
Jo 

a+pj^„i+i{ti)a^b{j,m)2cos[rt+{m){t - ti)] 
+CT-Pj,m-i(ti)CT+6(j, -m)2cos[ri_(TO)(t - ti)] 

-p„n ih b{j, m)e-'''^ (-^)(t~t,) 



(60) 



where h{j,m) = j(j + 1) — m(m + 1). Of course, the 
TCL version is obtained again by replacing pjm{ti) with 
Pjm{t) under the time integral. 

The procedure for the determination of the coherences 
and the populations of the central spin is analogous to 
the one of the previous case. Again, we find that the 
NZ and the TCL approach yield similar results for the 
coherences. An example is shown in Fig. [7] for the TCL 
approximation only. We see that the second order with 
the projection ([59|) leads to an excellent approximation 
for the coherences even for very long integration times. 



1 
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FIG. 7: (Color online) Coherence of the central spin. Top: 
Second-order TCL approximation for the projection ((59}. 
Bottom: Exact solution. Parameters: = 101 and a = 0.1. 



A further remarkable feature of the projection (|59p is 
that the corresponding NZ equation of second order ([50]) 
reproduces the exact dynamics of the populations. This 
can be demonstrated by solving the equations for the 
populations obtained from the master equation ([50]) . Al- 
ternatively, we can prove this statement directly with the 
help of the following argument. We first note that by use 
of Eq. ([5]) the von Neumann equation ([2]) leads to the 
exact equation 



from which we obtain an exact equation for the popula- 
tions of the central spin: 

|(±|Pp(t)|±> = ^* dt^{±\VmL{t^)p(t{)\±). (62) 

The crucial point is that for arbitrary p(t\) the following 
relation holds, 

(±|7'/:(i)/:(ii)p(ii)|±) = {±\vL{t)L{t{)Vp(t^)\±). 

(63) 

Hence, when calculating the diagonal elements of the 
density matrix we may insert at any time a factor of V in 
front of the density matrix p(t\). This fact is obviously 
related to the invariance of the subspaces spanned by the 
states |-)-) (g) |j, m) and |— ) ® |j, to-|- 1), a fact that has al- 
ready been used to solve the Schrodinger equation. For a 
formal proof one uses the expression l|59p to demonstrate 
that Eq. ([55]) is a direct consequence of the relation 



£t(ti)£t(i) (|±)(±| n,™) = ^ «.n,,„, (64) 



with certain system operators Ajm' ■ 

Using now the relation ((63|) in Eq. ([62]) we find 



d 
~dt 



±\Vp{t)\±)^ f dti{±\VC{t)C{ti)rp{h)\±). (65) 
Jo 



dtiVC{t)C{ti)p{ti 



(61) 



This equation is exact and coincides with the equation 
for the population that is obtained from the second-order 
NZ equation We conclude that the NZ equation 

of second order indeed yields the exact dynamics of the 
populations. 



IV. CONCLUSIONS 

We have investigated the performance of the corre- 
lated projection operator method in the case of a spin- 
bath model where a central spin couples to a spin bath 
via full Heisenberg interaction. The coupling constants 
have been chosen to be uniform such that one can con- 
struct an exact solution of the von Neumann equation. 
We have considered two different projection superopera- 
tors, one of which projects onto the eigenspaces of the 3- 
component J3 of the total bath angular momentum. The 
other projection operator projects onto the simultaneous 
eigenspaces of J3 and the square of the total bath an- 
gular momentum. The choice of these projections is mo- 
tivated by the fact that both Jg"* and are conserved 
under the dynamics given by the full system Hamilto- 
nian. The projection superoperators studied here thus 
fully exploit the known symmetries of the model, grant- 
ing the invariance of the expectation values of the con- 
served quantities. For both projectors we have solved 
the Nakajima-Zwanzig and the time-convolutionless mas- 
ter equation to second order in the interaction and have 
compared the results to the exact solution. 
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Both the NZ and the TCL master equation lead to a 
relaxation and deeoherence behavior that approximates 
the exact solution very well. Remarkably, also the posi- 
tion and the width of the partial revivals of the coherence 
and the populations could be reproduced. For all tested 
parameter sets, the TCL and NZ solutions were hardly 
distinguishable. In these cases the TCL approach is to 
be favored because of its much simpler structure. 

A remarkable property of the projection superopera- 
tor that projects onto the simultaneous eigenspaces of 
and J3 is that the corresponding second-order NZ equa- 
tion results in the exact dynamics for the populations. 
While the standard product state projection gives only a 
poor approximation, the above correlated projection re- 
produces the exact solution already in lowest order. This 
illustrates an important feature of the correlated projec- 
tion operator method. Namely, different choices for the 
projection V yield dynamic equations of different struc- 
ture for different sets of relevant variables. Therefore, 
the resulting expressions for the quantities of interest, 
e.g. the populations, may differ in all orders of the cou- 
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